Effects of boundary roughness on a Q-factor of whispering-gallery-mode lasing 

microdisk cavities 
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We perform numerical studies of the effect of sidewall imperfections on the resonant state broad- 
ening of the optical microdisk cavities for lasing applications. We demonstrate that even small 
edge roughness (< A/30) causes a drastic degradation of high-Q whispering gallery (WG) mode 
resonances reducing their Q-values by many orders of magnitude. At the same time, low-Q WG 
resonances are rather insensitive to the surface roughness. The results of numerical simulation ob- 
tained using the scattering matrix technique, are analyzed and explained in terms of wave reflection 
at a curved dielectric interface combined with the examination of Poincare surface of sections in the 
classical ray picture. 
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During recent years significant experimental efforts 
were put forward towards investigation of laser emis- 
sion of dielectric and polymeric low-threshold microdisk 
cavitiesiSiMifiiSiLSiP . The high efficiency of lasing oper- 
ation in such devices is related to the existence of natu- 
ral cavity resonances known as whispering gallery (WG) 
modes. The origin of these resonances can be envisioned 
in a ray optic picture, wherein light is trapped inside the 
cavity through total internal reflections on the cavity-air 
boundary. 

One of the most important characteristics of cavity 
resonances is their quality factor (Q-f actor) defined as 
Q = 2ir* (Stored energy) /(Energy lost per cycle). The 
Q-factor of a microdisk cavity is mostly governed by a 
radiative leakage through the curved interface due to 
diffraction. An estimation of the Q-fa,ctov in an ideal 
disk cavity of a typical diameter d ~ 10/im for a typical 
WG resonance gives Q ~ 10 13 (see below, Eq. ©). At 
the same time, experimental measured values reported 
so far are typically in the range of 10 3 ~ 10 4 or lower. 
A reduction of a Q-factor may be attributed to a variety 
of reasons including side wall geometrical imperfection, 
inhomogeneity of the height and diffraction index of the 
disk, effects of coupling to the substrate or pedestal and 
others. Several experimental studies point out side wall 
imperfections as the main factor affecting the Q-factor of 
the cavitjt£*£*£. An indirect indication of the importance 
of this factor in disc microcavities is provided by the ob- 
servation that typical Q-factors of spheroidal microcavi- 
ties are several orders of magnitude higher than those of 
microdisk of comparable dimension a 1 ' 10 . This is believed 
to be due to superior quality of the microsphere surfaces 
where boundary scattering may be limited by thermal 
fluctuations of the surface only. Therefore, the effect of 
surface roughness appears to be of crucial importance for 
the design, tailoring and optimization of Q-values of las- 
ing mikrodisk cavities. To the best of our knowledge, this 
effect has not been considered to date and warrants an 
investigation. 

In order to compute the resonant states of a cavity of 
an arbitrary shape we develop a new approach based on 
the scattering matrix technique. The scattering matrix 



technique is widely used in analyse of waveguides** as 
well as in quantum mechanical simulations 12 . This tech- 
nique was also used for an analysis of resonant cavities for 
geometries when the analytical solution was available 13 . 
Note that because the problem at hand requires a fine dis- 
cretization of the geometry, commonly used time-domain 
finite difference methods 1 ^ would be prohibitively expen- 
sive in terms of both computer power and memory. While 
a detailed description of the calculations will be given 
elsewhere, we present here the essence of the method. 

We consider a two-dimensional cavity with the refrac- 
tion index n surrounded by air. Because the majority 
of experiments are performed only with the lowest trans- 
verse mode occupied, we neglect the transverse (z-) de- 
pendence of the field and thus limit ourself to the two- 
dimensional Helmholtz equation. We divide our system 
in an outer and an inner regions. In the outer region the 
refraction index n is independent of the coordinate and 
the solution to the Helmholtz equation can be written in 
polar coordinates in the form 



-f-OQ 

*= {A q H^{kr)+B q H^\kj 



(1) 



* = E z (H z ) for TM (TE)-modcs, are the 

Hankel functions of the first and second kind of order 
q describing respectively incoming and outgoing waves, 
k = lu I c = 2tt/X. 

We define the scattering matrix S in a standard 
fashioniiii^ii 3 ., B — SA, where A, B are column vectors 
composed of the expansion coefficients A q , B q for incom- 
ing and outgoing states in Eq. Q). The matrix element 
Sqi q gives the probability amplitude of scattering from 
an incoming state q into an outgoing state q' . In order 
to apply the scattering matrix technique we divide the 
inner region into N narrow concentric rings. At each 
i-th boundary between the rings we introduce the scat- 
tering matrix S 1 that relates the states propagating (or 
decaying) towards the boundary, with those propagating 
(or decaying) away of the boundary. The matrices S 1 are 
derived using the requirement of continuity of the tangen- 
tial components for the E- and ii/-neld at the boundary 
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between the two dielectric media. Successively combin- 
ing the scattering matrixes for all the boundariesii*i£, 
S 1 ® . . . S N , we can relate the combined matrix to the 
scattering matrix S. 

To identify the resonant states of a resonant cav- 
ity we introduce the Wigner-Smith time-delay matrix 
Q = - (dSydk) S; 2 i 13 i 15 where the diagonal elements 
Q qq give a time delay experienced by the wave inci- 
dent in g-th channel and scattered into all other chan- 
nels. The Q- value of the cavity is Q = ujT£>(k), where 
To(k) is the total time delay averaged over all M incom- 
ing channela 2 i 13 i 16 , 
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exp(z# M ) = A M are the eigenvalues of the scattering ma- 
of the matrix S, det S 



trix S, = £ 



is the total phase of the determinant 
exp(i9). 
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FIG. 1: (a) Dependence of the quality factor Q = Q(X) of 
the circular disk for different surface roughness Ar indicated 
in the figure; the disk radius is R — 5/im, the refraction in- 
dex n — 1.8. The inset illustrates a cavity where the surface 
roughness Ar = 200nm and Al = 2irR/50 (the dotted line 
represents an ideal circular boundary, the shaded region cor- 
responds to the cavity) . p characterizes the average radius of 
local curvature due to boundary imperfections. TE modes of 
the cavity exhibit similar features and are not shown here. 

Figure^shows calculated Q- values of the disk resonant 
cavity for different surface roughnesses for TM modes in 
some representative wavelength interval. Note that an 
exact experimental shape of the cavity-surface interface 
is not available. We thus model the interface shape as 
a superposition of random Gaussian deviations from an 
ideal circle of radius R with a maximal amplitude Ar/2 
and a characteristic distance between the deviation max- 
ima AZ ~ 2-7T-R/50 (see illustration in inset to Fig. 2J. 

The solid curve in Fig. ^ corresponds to an ideal disk 
cavity without imperfections. Resonant states of an ideal 
disk (as well as the bound states of the corresponding 
closed resonator) are characterized by two numbers, q 



(see Eq. and m. The index m is a radial wave num- 
ber that is related to the number of nodes of the field 
components in the radial direction r. The angular wave 
number q can be related to the angle of incidence x m a 
classical ray picture^ 



q = nkR sinx- 



(3) 
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FIG. 2: Dependence Q on the surface roughness Ar for sev- 
eral representative resonances. (Each curve remains practi- 
cally unchanged for different realizations of surface rough- 
ness). The inset shows the dependence of local radius of 
roughness curvature p subject to Ar. Parameters of the cavity 
are the same as in Fig. 

The dependence of the averaged Q-values on the sur- 
face roughness Ar is summarized in Fig. [21 for several 
representative resonances. A common feature of all high- 
Q resonances is a dramatic decrease of their maximal Q- 
value that takes place for very small values of Ar < A/20. 
For example, a Q- value of at he resonant state TMs3,i 
drops from Q rs 10 13 for an ideal disk to Q « 10 3 for 
surface roughness of only Ar = 20nm. In contrast, low- 
Q resonances show a rather slow decrease in their Q- 
values over the range of variation of Ar. For example, 
for the same surface roughness Ar = 20 the Q-value of 
the resonant state TM^^ decreases only by a factor of 
1.5, dropping to Q k 200. 

In order to understand these features we combine a 
Poincarc surface of section (SoS) method with an anal- 
ysis of ray reflection at a curved dielectric interface 19 . 
The Q-value of the cavity can be related to the transmis- 
sion probability T of an electromagnetic wave incident 
on a curved interface of radius p by Q — 2nkpcosx/T— 
(this expression is valid for large angles of incidence \ 
when T <C 1). In turn, for knp 3> 1, the transmission 
probability readsi^ 



T = \Tp \ exp 



2 nkp 



3sin 2 (x) 



(cos 2 Xc ~ cos 2 x) 



3/2 



(4) 



where Tp is the classical Fresnel transmission coefficient 
for an electromagnetic wave incident on a flat surface, 
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Xc = arcsin(l/n) is an angle of total internal reflection. 
Figure illustrates that T decreases exponentially as the 
difference X ~ Xc grows. 

The inset to Fig. [21 depicts the Poincare SoS for two 
states with q — 56 and 83 shown in Fig. ^ where the 
initial angle of incidence xo of launched rays is related to 
the angular number q by Eq. J2J). The SoS demonstrates 
that initially regular dynamics of an ideal cavity trans- 
forms into a chaotic one even for a cavity with maximum 
roughness Ar < 20nm. AX^' 1 in Fig. ^indicates the es- 
timated increase in the transmission coefficient due to the 
broadening of the phase space, Axch, as extracted from 
the Poincare SoS for the state with q — 83. This cor- 
responds to the decrease of AQ ~ AT -1 w 10~ 2 . This 
value is much smaller that the actual calculated decrease 
of the Q-factor for the high-Q resonance TM^i. 

To explain the rapid degradation of high- Q resonances, 
we concentrate on another aspect of the wave dynamics. 
Namely, the imperfections at the surface boundary effec- 
tively introduce a local radius of surface curvature p that 
is distinct from the disk radius R (see illustration in Fig. 
^) . One may thus expect that with the presence of a local 
surface curvature, the total transmission coefficient will 
be determined by the averaged value of p rather than by 
the disk radius R. The dependence of p on surface rough- 
ness Ar for the present model of surface imperfections is 
shown in the inset to Fig. [21 Figure demonstrates 
that the reduction of the local radius of curvature from 
5/zm (ideal disk) to 1.7pm (Ar = 20nm) causes an in- 
crease of the transmission coefficient by AT cur as 10 s . 
This estimate, combined with the estimate based on the 
change of AT c h is fully consistent with the Q-factor de- 
crease shown in Figs. If 121 We thus conclude that the 
main mechanism responsible for the rapid degradation of 
high-Q resonances in non-ideal cavities is the enhanced 
radiative decay through the curved surface because the 
effective local radius (given by the surface roughness) is 
smaller that the disk radius R. 

For the case of low-Q resonances the change in the 
transmission coefficient due to enhanced radiative decay 
AT cur is of the same magnitude as the change AT c h due 
to the broadening of the phase space caused by the tran- 
sition to chaotic dynamics (as illustrated in Fig. [3] for 
the resonance TM^^). Therefore, both these factor play 
comparable roles in degradation of the low-Q WG reso- 
nances. 

It is worth mentioning that one often assumes that 



long-lived high-Q resonances in idealized cavities (e.g. 
in ideal disks, hexagons, etc.) are not important for 
potential application in optical communication or laser 
devicesi&l& because of their extremely narrow width. 
Our simulations demonstrate that it is not the case, be- 
cause in real structures the Q-values of these resonances 
becomes comparable to those of intermediate- Q reso- 
nances already for small or moderate surface roughness 
of Ar ~ 10 - 50 nm. 

To conclude, our results highlight the importance of 
surface roughness for the performance of microcavities 
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FIG. 3: Dependence T = T{\) for several radii of curvature 
p according to Eq. 0] Inset shows a Poincare SoS for the 
states q — 83 and q = 56 for the cavity with Ar = (straight 
lines of x = const) and Ar = 20nm. The number of bounces 
for a given angle of incidence xo is chosen in such a way that 
the total path of the ray does not exceed the one extracted 
from the numerically calculated Q- value for the corresponding 
resonance, L — ctd = Q/k. 



for laser applications, and provide estimations on surface 
roughness that might be instrumental for device design 
and fabrication. 
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